How to Track Protists in Three Dimensions 
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We present an apparatus optimized for tracking swimming microorganisms in the size range 
10—1000 fim, in three dimensions (3D), far from surfaces, and with negligible background convective 
fluid motion. CCD cameras attached to two long working distance microscopes synchronously 
image the sample from two perpendicular directions, with narrowband dark-field or bright-field 
illumination chosen to avoid triggering a phototactic response. The images from the two cameras 
can be combined to yield 3D tracks of the organism. Using additional, highly directional broad- 
spectrum illumination with millisecond timing control the phototactic trajectories in 3D of organisms 
ranging from Chlamydomonas to Volvox can be studied in detail. Surface-mediated hydrodynamic 
interactions can also be investigated without convective interference. Minimal modifications to the 
apparatus allow for studies of chemotaxis and other taxes. 

PACS numbers: 87.17. Jj,87.18.Ed,47.63.Gd 
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I. INTRODUCTION 

Protists, the grouping of eukaryotic microorganisms 
that encompasses such diverse entities as flagellated and 
ciliated protozoa |l| (e.g., Euglena, Paramecium) and 
motile green alga [2| (e.g., Chlamydomonas, Volvox - 
shown in Fig. [TJ, constitute an important class of or- 
ganisms in the study of evolutionary biology, biological 
physics, and, recently, biological fluid dynamics 0,0, @. 
Many flagellated protists display swimming behavior that 
is inherently three-dimensional (3D). A number of im- 
portant questions in biology and physics are associated 
with how the motion of such organisms is related to their 
body plan and to external stimuli such as light 0, 0] , dis- 
solved molecular species, gravity, temperature, bound- 
aries [8], and electromagnetic fields. It is thus desir- 
able to track their position and orientation in 3D 0, [l(| 
with high spatio-temporal resolution and, unless desired, 
free from systematic bias introduced by external stim- 
uli, background fluid motion, and hydrodynamic surface 
effects [ll 03. 

The first apparatus able to track microorganisms in 
3D was designed for bacteria [l3| and utilized an ana- 
logue feedback loop that moved the microscope stage to 
keep a bacterium centered in the field of view. Larger 
microorganisms such as protists require larger sample 
chambers, reaching millimeters or even centimeters in 
depth to avoid boundary-induced hydrodynamic effects. 
In this regime, methods based on a moving microscope 
stage are not suitable, as they induce uncontrolled back- 
ground fluid motion in the sample chamber. Similar con- 
siderations enter the study of the millimetric nematode 
C. elegans crawling on the surface of agar, for which a 
moving substrate introduces unwanted mechanical stim- 
uli. Instead, the camera itself can by moved by mo- 
tors controlled by an algorithm that dynamically centers 
the worm in the field of view [l4j |. Several apparati for 
tracking the 3D motion of microscopic particles without 
a moving stage have emerged in recent years, yet none 
is ideal for studying the swimming of protists. Meth- 



ods based on controlled defocussing of particles [H, EBl , 
placing a cylindrical lens in the imaging optics of the 
microscope [l3], or measuring the deflections of a laser 
beam that is focused close to a particle [ID, [l!| all suffer 
from a small tracking range along the optical axis. Ob- 
serving a sample which is illuminated from the side with 
a continuous gradient of color in order to color-code the 
third dimension [2(3, [H| suffers from low spatial resolu- 
tion in that dimension, and may provide a photostimulus 
to protists. Tracking objects with a confocal microscope 
is only possible when the objects move at very low speeds 
[22l [23 | . Digital in-line holography may also be used for 
3D particle tracking, yet even vibrations with amplitudes 
< 1 (tzm of components along the optical path lead to a 
time-varying background in the hologram that can sig- 
nificantly degrade the signal of the moving object. 

The difficulties mentioned above can be overcome by 
using more than one camera to observe synchronously 
the sample chamber from different angles and then com- 
bining the images to yield 3D tracks of particles. Yet, 
existing implementations of this technique can not be 
easily adapted to track microscopic objects [HI HE]; t ne > 
either have a controlled but undesirable shear flow [2i 
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FIG. 1: Two protists whose swimming motion is of interest in 
this work. (A) Chlamydomonas reinhardtii (scale bar 5 //m) 
and (B) Volvox carteri (scale bar 200 pirn). 



have not integrated any control of thermal convection 
[27| . or have been implemented with a temperature gra- 
dient across the sample (~ 3 K/cm) to eliminate thermal 
convection 28]. The latter method is an important ad- 
vance, but may introduce an unwanted behavioral stim- 
ulus to protists. 

Here we present an apparatus that uses two identical 
imaging assemblies at right angles to each other that can 
be used in dark- and bright-field illumination, combined 
with systems to deliver photo-stimuli, to control tem- 
perature inside the sample, and to eliminate background 
fluid motion in large sample chambers (up to 2.5x2.5x5 
cm) filled with aqueous growth medium, by homogeniz- 
ing the temperature inside the sample chamber to mil- 
likelvin precision. The apparatus uses almost entirely 
off-the-shelf components and requires minimal expertise 
in optics. Online supporting material includes software 
to control the hardware and to perform 3D tracking. We 
discuss the resolution and limitations of the apparatus, 
present swimming trajectories of the protists Chlamy- 
domonas and Volvox in 3D, and illustrate with dual- view 
particle imaging velocimetry (PIV) the flow field that 
Volvox generates near a surface. 



II. EXPERIMENTAL APPARATUS 

The 3D tracking system, shown schematically in Fig. 
[H is based on a flexible but powerful imaging system, 
phototactic stimulus lights, and equipment to control and 
homogenize the temperature inside the sample chamber. 
These three elements are now explained in detail. 

The imaging system is comprised of two identical as- 
semblies that are mounted at right angles on a vibra- 
tion isolation table (Science Desk, with 900 x 1200 x 60 
mm breadboard, Thorlabs, Ely, UK). A monochrome 
Fire Wire CCD camera (Pike F145B, Allied Vision Tech- 
nologies, Stadtroda, Germany; 1388 x 1038 pixels, each 
6.45 x 6.45 /mi, maximal frame rate of 30 fps, and support 
for external triggering) was attached to each of the two 
microscopes (InfiniVar CFM-2/S, Infinity Photo-Optical, 
Boulder, CO). These are continuously focusable with a 
working distance between 18 mm and oo, yielding a max- 
imum magnification of x 9 at the smallest available work- 
ing distance. To allow a variable working distance, the 
camera/microscope assemblies were mounted on sliding 
rails (PRL-12, Newport Corp., Irvine, CA) via standard 
post/post support hardware. The horizontal rail was at- 
tached directly to the breadboard while the vertical one 
was attached to a movable and lockable rail carrier on a 
large optical rail (X95, Newport Corp.) mounted verti- 
cally to the optical breadboard. The outer chamber (see 
Fig. [2j and details below) limits the smallest working 
distance to ~ 60 mm, yielding a magnification of x 1 . 
While sufficient to image protists of size ~ 10 /jm with 
dark-field illumination, such organisms are better visual- 
ized with an additional x2 magnifier lens (2xDL, Infinity 
Photo-Optical). This increases the working distance and 




FIG. 2: Schematic drawing of apparatus. The outer chamber 
(light blue, lid and flange in grey) contains a water bath, fed 
through inlets and outlets (purple), and mixed with a mag- 
netic stir bar (white) driven by a motor external to the tank 
(not shown). The sample chamber (white) is suspended by 
stainless steel holders (light green), and illuminated by annu- 
lar LED arrays (red) . Microorganisms are visualized with two 
long-working-distance microscopes (dark gray) equipped with 
CCD cameras (dark green). Phototactic stimulus is provided 
by two LED and lens assemblies (yellow), and controlled by 
shutters (orange). 



the depth of field, as discussed below. 

The imaging system can be used with dark- or bright- 
field illumination: bright-field can be advantageous for 
organisms larger than ~ 100 /im as it captures more 
details of the organism (e.g., the body axis); dark-field 
can be desirable when the organism is so small that only 
center-tracking is possible, as under these conditions it 
yields a better signal-to-noise ratio. The flexibility to 
use both illumination techniques is obtained by using 
an annular LED array (LFR-100-R, CCS Inc., Kyoto, 
Japan) as the (unpolarizcd) light source for each micro- 
scope. Dark-field illumination is achieved when the field 
of view of the microscope only includes the dark region 
in the center of the LED annulus. Bright-field illumina- 
tion is achieved by inserting a diffuser plate (bulk frosted 
acrylic cut to size, RS components, UK) in between the 
LED array and the outer chamber, as far from the LED 
array as possible. We chose the color of the LED array to 
be narrowband red (655 nm, 21 nm bandwidth) as it has 
been shown that this color does not trigger a phototactic 
response in motile green algae [29| . 

For phototaxis studies, two opposing light sources are 
required in order to observe reproducible light-induced 
U-turns. The photo-stimulus lights were two broad spec- 
trum cool white Luxeon LEDs (MWLED, Thorlabs), col- 
limated with a short focal length lens (/ = 35 mm, 
LA1027, Thorlabs), and mounted to a shutter with mil- 
lisecond precision (SH05, Thorlabs). The shutter/LED 
assemblies are mounted via standard cage hardware 
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(Thorlabs) to posts on rotating stages (RP01, Thorlabs), 
allowing the direction of the light to be controlled. The 
stimulus direction was typically chosen to be horizontal 
in order to avoid an additional gravitational stimulus. 
This choice forces the common axis of the two cameras 
to be horizontal, along the stimulus direction. The beam 
diameter was controlled with an iris (SM1D12, Thorlabs) 
to illuminate only those faces of the sample chamber 
perpendicular to the common axis of the two cameras, 
thereby avoiding reflections, and the resulting unclear 
stimuli, from the other four faces. 

The two CCD cameras and shutters for the photo- 
stimulus lights were controlled with labview including 
the add-on toolbox ni-imaq for 1394-ieee cameras (Na- 
tional Instruments, Austin, TX), allowing precise syn- 
chronization of image acquisition. A labview program 
is part of the supporting online material. 

In order to combine the images from both cameras to 
yield 3D swimming tracks, it is crucial that both micro- 
scopes operate at the same magnification (i.e. the same 
working distance). This is easily achieved to sufficient 
precision by replacing the sample chamber with a tilted 
calibration ruler that is observable through both cam- 
eras, and adjusting the working distance until the field 
of view of both cameras has the same physical size. After 
this calibration step, the microscopes need to be aligned 
along their common axis such that the field of view of 
both cameras contains the same section of the common 
axis. This axis alignment greatly improves the ability to 
reconstruct 3D tracks, as explained in Section IIIII 

To obtain swimming trajectories that are not influ- 
enced by hydrodynamic surface effects or background 
flows, it is desirable to have a sample chamber that is 
as large as possible, while maintaining the fluid within it 
perfectly still. A stationary fluid can only be obtained 
if the temperature in the chamber and of the chamber 
walls is very homogeneous, thereby eliminating thermal 
convection caused by heating from the two LED arrays 
(each LED array consumes ~ 3.6 W). For a closed cham- 
ber, the critical Rayleigh number above which thermal 
convection starts to occur is Ra c = agL 3 AT/i/K ~ 1708, 
where a is the thermal expansion coefficient of the fluid, 
g is the gravitational acceleration, n is the thermal dif- 
fusivity, v is the kinematic viscosity, and L is the length 
scale across which there is a temperature difference AT 
[3fj| . While the precise value of Ra c depends on the geom- 
etry of forcing and on boundary conditions, the scale of 
temperature differences involved for water (a = 2x 10~ 4 
K _1 , k = 1.4 x 10~ 3 cm 2 /s, v = 0.01 cm 2 /s) is roughly 
100 mK for a one centimeter length. The largest cham- 
bers that have previously been temperature-homogenized 
below the thermal convection threshold, under compara- 
ble conditions to those presented here, have L ~ 1 cm 
in sedimentation studies [H], H3 ■ The system presented 
here eliminates thermal convection in chambers as large 
as 2.5 x 2.5 x 5 cm, implying temperature differences be- 
tween faces of the chamber to be below ~ 8 mK. 

Recognizing that in Stokes flow the effects of bound- 
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FIG. 3: Decay of convective motion in the inner sample cham- 
ber. Curves show the root-mean-squared velocity « rms within 
the chamber, obtained by PIV, as a function of time for three 
different chamber dimensions: blue - 2.5 x 2.5 x 2.5 cm, ma- 
genta - 2.0 x 2.0 x 2.0 cm, green - 1.5 x 1.5 x 1.5 cm. Inset 
shows the time constant for each decay as a function of cham- 
ber size, consistent with the L 2 diffusive scaling. 



aries at a distance h from compact objects acted upon by 
gravity fall off as h , and given that swimming trajecto- 
ries can easily sample a vertical scale that is 5 times the 
organism diameter, the sample chamber should be > 20 
times the organism diameter for surface effects to remain 
below the 5% level [IJ G3- A chamber of this size thus 
allows protists as large as 1000 /jm to be studied with 
negligible hydrodynamic surface effects. 

Before explaining how the temperature in the sample 
chamber is controlled and homogenized, it is necessary to 
give details of the outer and sample chambers. The outer 
chamber has dimensions 12 x 12 x 10 cm, two inlets and 
two outlets (as shown in Fig. \2§ and is made from 2.75 
mm thick borosilicate glass (custom made by Fine Glass 
Finishers Ltd., Great Chesterford, UK). The flange and 
lid of the outer chamber were custom made out of PVC 
with a CNC machine. The three types of sample chamber 
used were (i) custom cuvettes made from standard micro- 
scope slides cut with a tungsten glass scriber (LAC-450- 
A, Fisher Scientific, UK) and glued together with UV- 
curing optical glue (NOA68, Norland Products, Cran- 
bury, NJ), cured in a UV Chamber (ELC-500, Electro- 
Lite Corporation, Bethel, CT), (ii) non-standard com- 
mercial glass cuvettes (VitroCom, Mountain Lakes, NJ), 
and (iii) standard glass cuvettes. The sample chamber is 
held rigidly in the center of the outer chamber by four 
thin stainless steel rods (1.6 mm diameter). Each rod was 
inserted into a corresponding tight-fitting hole in the lid 
and fixed with a set screw, allowing easy modification of 
the holding arrangements for the different chambers. 

It is well-known that the flagella of protists such as 
Chlamydomonas have a strong tendency to stick to glass 
surfaces [33|. As cells swim around the chamber they 
inevitably collide with the chamber walls. To avoid any 
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problem with sticking, the glass is coated with PDMS 
(Sylgard 184, Dow-Corning, Belgium) and etched for 3 
minutes in a plasma cleaner (Femto, Diener Electronic, 
Nagold, Germany), following a published protocol [34j . 

The temperature in the sample chamber is controlled 
by cycling filtered water through the outer chamber from 
a large tank that contains a submersed thermostatic 
heater (Jager Eheim, Deizisau, Germany) and pump 
(Eheim, Deizisau, Germany). In order to homogenize 
the temperature in the sample chamber, the pump is 
switched off, and strong rare earth magnets (EP200, e- 
Magnets UK) are spun at ~ 300 rpm by a sturdy mo- 
tor (178-5112, RS components, UK) on the outside of 
the outer chamber, thereby moving a Teflon-coated 5 cm 
magnetic stir bar inside the outer chamber at the same 
speed. This stirring evens out the temperature within 
the outer chamber and, if the stir bar is spun at the ap- 
propriate speed and place in the outer chamber, moves 
the water past the faces of the sample chamber with- 
out setting up recirculating vortices on the faces. By 
injecting inexpensive tracer particles (size <75 fj,m, Pli- 
olite VTAC-L, Eliokem, Villejust, France) into the outer 
chamber in order to visualize the flow across the faces 
of the sample chamber, we found that the arrangement 
drawn in Fig. [5] can homogenize the temperature inside 
the sample chamber below the threshold Ra c . 

The flows within the sample chamber were quantified 
with commercial PIV software (FlowManager, Dantec 
Dynamics, Skovlunde, Denmark). As a metric for the 
extent of flows throughout the chamber, we report the 
r.m.s. velocity v ims obtained by uniform averaging over 
the field of view of one camera. Figure [3] shows that 
v Tms decays exponentially as a function of time after the 
stirring of the water bath has begun, with a time con- 
stant t(L) that depends on the smaller chamber dimen- 
sion L. It is straightforward to see when the tempera- 
ture difference falls below the threshold for convection, as 
tracer particles that are used to visualize the convective 
flow in the sample chamber will suddenly begin to fall 
out of the fluid at approximately their Stokes sedimen- 
tation speed, forming a sedimentation front that propa- 
gates downward. For 10 ^m latex beads (C37259, Invit- 
rogen, Carlsbad, CA), the Stokes sedimentation speed is 
~ 3 /im/s. This serves as a lower bound for v rms in Fig. 
[3j We expect t(L) to arise from viscous dissipation, and 
thus to scale as t(L) ~ L 2 jv. In water this yields times 
on the order of a few to ten minutes for L in the range 
1 — 2.5 cm, consistent with the data in Fig. [3J Each curve 
conforms well to a single exponential decay, and the fit- 
ted times t{L) obey well the expected quadratic scaling 
as shown in the inset to Fig. [3l with a value not far from 
that expected from the kinematic viscosity of water. 



III. TRACKING SOFTWARE 

The 3D tracking was done by analyzing the image 
sequences from each camera separately, giving a set of 



two two-dimensional (2D) tracks, and combining suitable 
tracks from these two sets to yield 3D trajectories. 

Modified matlab (MathWorks, Natick, MA) versions 
of freely available [351 ] particle tracking routines written 
by J.C. Crocker and D.G. Grier were used for the 2D 
tracking, allowing many organisms to be tracked at the 
same time. To track organisms that are so small that 
they appear circular and without internal structure in 
dark-field images (e.g., Chlamydomonas) , no modifica- 
tions to the original versions of the code need to be made. 
To track extended objects (e.g., Volvox) the routines that 
identify the center of the object need to be modified in a 
way that depends on the shape and structure of the ob- 
ject. For bright-field images of the spherical Volvox, all 
internal structure of Volvox was removed by histogram 
equalization, followed by spatial band pass filtering. The 
resulting image was then convolved with a binary disk- 
shaped kernel, yielding an image in which the centroids 
of peaks correspond to the Volvox centers in the original 
image. A Volvox colony carries daughter colonies inside 
it (see Fig. [T]), which are fixed in the posterior hemi- 
sphere and therefore act as convenient markers of the 
body axis. The axis of Volvox can thus be determined by 
finding the vector between the geometric center of Volvox 
and the center of brightness of the daughter colonies in- 
side Volvox for both directions and then combining these 
two vectors to obtain a 3D axis. The modified code also 
allows additional information for each Volvox to be gath- 
ered, such as the orientation of the body axis. 

To identify two 2D tracks that are suitable for synthe- 
sizing into a 3D track, the two sets of 2D tracks were com- 
pared along the common axis of the field of view of the 
two cameras. Consider one of the possible combinations 
of two 2D tracks. Even if these two tracks are projec- 
tions of positions of a single organism, the tracks usually 
do not completely overlap in time, because during the 
course of a long track the signal from the tracked object 
may fall below the tracking threshold so that the object 
'drops out' of the tracking data [36]. This means that 
only the time-overlapping sections of each track can be 
compared. Because of the precise alignment of the field 
of view of both cameras (see Sec. QI]), a decision upon 
whether the two 2D tracks are from the same organism 
can be made by finding the r.m.s. difference between the 
position-coordinate along the common axis. The two 2D 
tracks for which this value is minimal (and below a cer- 
tain threshold) arc then synthesized into a 3D trajectory. 
Code that can perform all the operations described above 
is part of the supporting online material. 



IV. PERFORMANCE 

The tracking precision of the apparatus has been de- 
termined by the standard method [3?], [38| of observing 
fluctuations in the tracked position of particles that are 
fixed between two cover slips. The precision was tested at 
the minimum working distance the apparatus allows (60 
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mm, corresponding to x 1 magnification) , for two differ- 
ent types of objects. For monodisperse 10 /im latex beads 
imaged in dark-field illumination, the uncertainty in the 
position was found to be < 1.5 fim, if the x2 magnifier 
lens is mounted to the microscope. For objects in the size 
range 425-500 /im (Volvox fixed with iodine), in bright- 
field illumination, the uncertainty was found to be < 1.3 
/mxi without the magnifier lens. 



The performance of the apparatus is determined not 
only by the uncertainty in spatial position, but also by 
the overall volume in which 3D tracking can be per- 
formed. This volume is set by the depth of field of the 
microscope. For the purpose of simply tracking the cen- 
ter of a spherical object, the object can be out of focus 
as long as the signal in the image intensity profile is suf- 
ficiently large. Therefore the "trackable depth" (TD), 
which we define as the depth in which the signal/noise 
< 4, is a more suitable measure for the trackable volume 
than the depth of field. The TD is strongly dependent on 
the object size and on the signal (dark-field illumination 
gives a larger signal). We measured the TD for 10 /xm 
beads in dark-field to be 5.8 mm (11.9 mm) at a work- 
ing distance of 60 mm (100 mm). Imaged in bright-field, 
fixed Volvox of size 425-500 /jm have TD = 18.3 mm at 
the minimum working distance of 60 mm. 



A limitation of the tracking software presented here is 
that the concentration of organisms in the sample cham- 
ber should not be so large that swimmers overlap fre- 
quently in the 2D images from each camera. As the 
tracking software can not distinguish overlapping objects, 
a high concentration of swimmers would result in very 
short 2D tracks, and therefore less accurate synthesizing 
of 2D tracks to a 3D track. An alternative method for 
obtaining 3D tracks from the images of more than one 
camera is to determine the 3D position of every parti- 
cle at every time point. This approach is often taken in 
3D Lagrangian particle tracking [U and can han- 
dle larger concentrations of trackable objects, but is usu- 
ally implemented with at least three cameras in order to 
reduce frequent ambiguities in the 3D particle identifi- 
cation, and requires an elaborate calibration. This 3D 
tracking method could be implemented with the appara- 
tus described here, if one of the stimulus lights is replaced 
by a third microscope-camera-assembly. 



Another limitation is the size of the outer chamber, 
which limits the magnification to values that are not 
sufficient for tracking small bacteria (e.g., E. colt), even 
though the microscope has a maximum magnification of 
xl8 (including the x2 magnifier lens). For trac king bac- 
teria sample chambers of size 2x2x2 mm (401 may 
be used, for which the temperature-homogenizing outer 
chamber is not needed. 
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FIG. 4: Reconstruction of a swimming trajectory of Chlamy- 
domonas reinhardtii. Gravity is along negative ^-direction. 
(A) 3D trajectory, color-coded to indicate local speed. Red 
dot indicates start of the trajectory. (B)-(D) Three compo- 
nents of position versus time. In (B) are the two traces x(t) 
from the two cameras. 



V. APPLICATION OF APPARATUS TO 
TRACK GHLAMYDOMONAS AND VOLVOX 

The apparatus was tested on two low- Reynolds number 
swimmers of very different size: Chlamydomonas rein- 
hardtii (diameter of ~ 10 fim) and Volvox barberi (diam- 
eter of ~ 600 /im). Both species we re g rown axenically 
in Standard Volvox Medium (SVM) [4JJ with sterile air 
bubbling, in diurnal growth chambers (Binder KBW400, 
Tuttlingen, Germany) set to a daily cycle of 16 h in cool 
white light (~ 4000 lux) at 28° C and 8 h in the dark at 
26° C. Sample chambers were filled with SVM and were 
of size 2.5 x 2.5 x 5 cm for Volvox, and 1x1x4 cm 
(standard cuvette) for Chlamydomonas. 

The biflagcllatcd Chlamydomonas beats its flagella at 
~ 40 Hz, primarily in the manner of the breast stroke. 
Its most familiar swimming trajectory is helical, with a 
radius of 20 /im and a speed on the order of 50 //m/s. 
The cell has an "eye-spot" that serves as a photosensor, 
and the changing illumination levels of the eye-spot lead 



FIG. 5: A phototactic turn of Volvox barberi. (A) The posi- 
tional and orientational measurements are illustrated by vec- 
tors indicating the body axis, and swimming-speed-dependent 
coloration of the track. To initiate the 180° change swimming 
direction, the light initially was from the right along the ir-axis 
and then changed to come from the left. Gravity is directed 
along the negative z-direction. Sphere represents initial posi- 
tion along the track. (B) Evolution of the body axis during 
the phototactic turn is described in terms of two angles 9 and 




to transient changes in flagellar beat dynamics in such a 
manner that the cell can turn toward the light. 

Figure UJA. shows a 45 s trajectory, obtained without 
phototactic stimulus, during which the cell explored a 
volume less than 1 mm 3 . The manner in which the two 
2D trajectories from the cameras are synthesized into a 
3D trajectory is indicated in panels (B)-(D) of the figure. 
The x-axis is common to the two cameras, and the over- 
lap between the two is clear in Fig. [4)3 . The very slight 
graded mismatch between the two ^-component curves 
reflects a slight misalignment of the cameras, and is eas- 
ily removed by remapping the pixel coordinates. 

Volvox barberi typically has 10, 000 — 50, 000 biflagel- 
lated somatic cells rigidly embedded at the surface of a 
transparent extracellular matrix. These beat at a typi- 
cal frequency of 20 Hz, primarily from the anterior pole 
to the posterior pole, with a slight tilt of the beat plane 
which leads to the characteristic spinning motion as it 
swims at speeds up to 800 /im/s [42j . Just as in Chlamy- 
domonas, each somatic cell has an eye-spot that modu- 
lates the beating of its two flagella, allowing the whole 
colony to perform phototaxis. A 3D track of the Volvox 
center and body axis during a phototactic turn is shown 
in Fig. Determining the body axis of Volvox by us- 
ing the position of the daughter colonies, as explained in 
Sec. IIII1 leads to the time series of the angles and <j> 
in Fig. 03. This is slightly noisy because the daughter 
colonies are not distributed evenly. The track also shows 



FIG. 6: Volvox carteri swimming near a surface. Flow fields 
from particle imaging velocimetry of a colony swimming up- 
ward against a horizontal cover slip that is glued into the 
sample chamber (of size 2.5 x 2.5 x 5 cm), as seen from the 
side (A), and the top (B). Images were taken at x2 magnifi- 
cation. Scale bar is 200 /im. 



an interesting balance between the bottom-heaviness of 
Volvox (due to the clustering of daughters in the poste- 
rior hemisphere), which tends to align the axis with the 
^-direction, and the phototactic tendency to align the 
colonial axis with the direction of the light (the x-axis). 

In addition to allowing a controlled systematic distur- 
bance to the behaviour of microorganisms in the form 
of light, this apparatus is also suitable to study the in- 
teractions of microorganisms with a surface, as the sam- 
ple chambers can be made so large and so temperature- 
homogenized that the effect from other surfaces and ther- 
mal convection can be neglected. One interesting ef- 
fect, discussed in more detail elsewhere occurs when 
colonies swim upward to a horizontal surface and rotate 
in place. As the colonies are denser than the water in 
which they swim, and thus have a net external force act- 
ing on them, the far-field flow around them is described 
by a stokeslet pointing downwards. In the neighborhood 
of a no-slip surface, the stokeslet induces a set of im- 
age singularities which together produce a characteristic 
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FIG. 7: Tracking of spinning V. carteri. Shown are the x 
(red) and y (blue ) coordinates of a colony spinning near an 
upper surface with its axis vertical. Recorded at 10 fps in 
convection controlled chamber. The oscillations represent a 
small periodic wobble in the colony centroid. 



lobed flow field [43|. This is readily demonstrated by 
glueing a horizontal surface (a microscope coverslip) into 
the sample chamber and performing dual-view PIV. The 
results of this are shown in Fig. HJ'VjB, illustrating the 
fluid flow vector fields viewed from both the side and the 
top. From above, we see vectors oriented predominantly 
inward. This inward flow leads to complex behavior of 
nearby colonies [8|. 

When viewed from above as in Fig. [6)3 this setup 
also provides a means to monitor the rotational dynam- 
ics of colonies in great detail, providing accurate mea- 
surements of the mean rotational frequency, the noise 
in rotational motion, and lateral drifts. As mentioned 
earlier, the bottom-heaviness of the colonies keeps the 
colonial axis oriented vertically, allowing the daughters 
to serve as convenient markers to track rotation. Pre- 
cise determination of time series of rotation can then be 
achieved by determining the correlation between succes- 
sive images and a reference image, adjusted for centroid 
drift. The centroid dynamics itself serves as a sensitive 



measure of colony asymmetries, such as mismatch be- 
tween the colonial axis and the axis defined by the center 
of buoyancy and the geometric center. Figure [7| shows 
the two components of the centroid position for a colony 
rotating against an upper surface, showing a clear peri- 
odic wobble at a frequency of ~ 0.5 Hz. The systematic 
drifts in position, which may be due to the swimming 
dynamics themselves or residual convective currents, are 
in any event below 2 /im/s. 



VI. CONCLUSIONS 

We presented an apparatus that can track swimming 
microorganisms in the size range 10 — 1000 /mi in 3D, 
without the influence of systematic bias due to beav- 
ioral stimuli, hydrodynamic interactions with surfaces, 
and convective background flows. As the apparatus can 
eliminate these biases, it can also be used to study the 
influence of each of them. The simplicity of the appa- 
ratus compared to other 3D tracking systems, and the 
software that is part of the supporting online material, 
make this system easily reproducible. 

The heart of the apparatus can also be used as the 
basis for studies of more complex phenomena. For in- 
stance, the entire device can be mounted on a tiltable 
platform in order to examine the effects of varying di- 
rection of gravity with respect to the phototactic axis. 
Likewise, a rotatable chamber can be substituted for the 
ones discussed here in order to examine the effects of fluid 



vorticity on phototactic swimming 44 1 . 
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